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METHODS FOR DETERMINING OPTIMAL DOSES FOR A FIRST 
INVESTIGATIONAL ADMINISTRATION IN HUMANS 



[001] CROSS-REFERENCE TO RELATED APPLICATIONS 
[002] The present application claims the benefit 
of United States Provisional Application No. 60/468,678, 
filed May 7, 2003, the entire disclosure whereof is 
incorporated herein by reference. 

[003] TECHNICAL FIELD OF THE INVENTION 
[004] The present invention relates to a method 
for determining an optimal dose range of a compound for a 
first investigational administration of that compound in 
a human. The method involves stochastic simulation of 
allometric parameters to produce outcome distributions 
that provide a basis for selecting an optimal dose range. 
The present invention also relates to optimal doses of 
compounds for use in the first investigational 
administration of such compounds to humans. 

[005] BACKGROUND OF THE INVENTION 

[006] A critical step in drug development is the 
transition from studies conducted on animals to in vivo 
testing in humans. The first administration of a 
compound under investigation to a human, as part of a 
first-in-human study (hereinafter " FIH study"), has to be 
conservative in dosage amount for safety reasons. Yet, 



the fiscal imperatives of drug development argue for an 
aggressive approach for testing in humans so that the 
first observable pharmacologic effect and the maximum 
tolerated dose levels are rapidly identified. This is 
especially true for compounds that will be given to a 
target human population during Phase I (e.g., anticancer 
agents, AIDS drugs, etc.). 

[007] Data obtained in the pharmacological and 
toxicological studies conducted on animals are a vital 
link in the prediction of pharmacokinetic /pharmacodynamic 
parameters in humans. The extrapolation of animal data 
to humans, known as interspecies scaling, is useful in 
translating dosage regimens from animals to humans. 
Interspecies scaling is based on the assumption of 
anatomical, physiological, and biochemical similarities 
across animal species. One common tool to perform 
interspecies scaling is allometric scaling. 

[008] Allometric scaling relates a 
pharmacokinetic parameter, e.g., clearance, to total body 
weight (BW) as follows: 

Y = a(BW) p (1) ; 

wherein Y is a pharmacokinetic parameter and BW is the 
body weight of the animal studied. cx and (3 are 
allometric parameters. Typically, allometric scaling is 
applied retrospectively to pre-clinical interspecies 
animal data in order to extract the allometric 
parameters. Then, the allometric parameters so extracted 
are applied prospectively to compute pharmacokinetic and 
pharmacodynamic responses in humans. See, e.g., "The 
Pharmacokinetic Principles Behind Scaling From 
Preclinical Results To Phase I Protocols," Mahmood, I., 
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and Balian, J. D., Cli . Pharmco . Kinet . , 1999 Jan: 36(1), 
incorporated herein by reference; "First-Time-in-Human 
Dose Selection: Allometric Thoughts and Perspectives, " 
Boxenbaum, H. and DiLea, C, J. Clin. Pharmacol., 1995; 
35:957-966, incorporated herein by reference; "Man Versus 
Beast: Pharmacokinetic Scaling In Mammals," Mordenti, J., 
Journal of Pharmaceutical Sciences, Vol. 75, No. 11, 
November 1996, incorporated herein by reference; 
"Interspecies Scaling And Comparisons In Drug Development 
And Toxicokinetics," Ings, R. M. J., Xenobiotica, 1990, 
Vol. 20, No. 11, pp. 1201-1231, incorporated herein by 
reference; "Computer-assisted drug development (CADD) an 
emerging technology for designing first- time-in-man and 
proof -of -concept studies from preclinical experiments," 
Gomeni, R. et al . , European Journal of Pharmaceutical 
Sciences, vol 13 (2001), pp. 261-270, incorporated herein 
by reference; "Interspecies Scaling, Allometry, 
Physiological Time, And The Ground Plan Of 
Pharmacokinetics," Boxenbaum, H. , Journal of 
Pharmacokinetics and Biopharmaceutics , Vol. 10, No. 2, 
pp. 201-227, 1982, incorporated herein by reference. 

[009] Typically, allometric scaling methods are 
applied deterministically to predict human 
pharmacokinetic parameters; i.e., the relationship 
between a pharmacokinetic parameter, such as clearance, 
has been treated as fixed and not subject to any 
uncertainty despite the known limitations of using animal 
data to predict human data. Consequently, any 
uncertainty in the allometric parameters determined from 
that relationship will lead to uncertainty in the 
extrapolation of animal data to predict human data. 

[0010] Thus, there is a need for a method that 
addresses the uncertainty in the estimates of the 
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allometric parameters. There is also a need for a 
method for determining optimal dose ranges for 
administering a compound under investigation for the 
first time in a human. 

[0011] SUMMARY OF THE INVENTION 

[0012] The present invention provides a method 
for determining the distribution of at least one outcome 
resulting from the administration of a compound to at 
least one human, comprising the steps of: 

(a) using non-linear mixed effects modeling to 
determine a value and a standard error for each of a 
plurality of allometric parameters ; 

(b) inputting said value and said standard error of 
at least one of said plurality of allometric parameters 
into a stochastic pharmacokinetic model, wherein said 
allometric parameter is designated as a random variable; 
and 

(c) using said stochastic pharmacokinetic model to 
computationally simulate administration of said compound 
to said human to produce said distribution of said 
outcome . 

[0013] The present invention also provides 
optimal doses of compounds under clinical evaluation in 
the first administration of such compounds to humans. 

[0014] DESCRIPTION OF THE FIGURES 
[0015] FIGURE 1 depicts the population 

pharmacokinetics as described by non-linear mixed effects 

allometric modeling. 

[0016] FIGURE 2 depicts the correlation between 

the predicted and the observed VX-7 02 concentration. 
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[0017] FIGURE 3 depicts the stochastic 
pharmacokinetic model used to computationally simulate 
the administration of VX-702 to human subjects. 

[0018] FIGURE 4 depicts the distribution of the 
predicted AUC for various dose levels. 

[0019] FIGURE 5 depicts the predicted 
distribution of average concentration for various dose 
levels . 

[0020] FIGURE 6 depicts the probability of the 
predicted median concentration for various dose levels. 

[0021] DETAILED DESCRIPTION OF THE INVENTION 
[0022] The present invention provides a method 
for determining the distribution of at least one outcome 
resulting from the administration of a compound to at 
least one human, comprising the steps of: 

(a) using non-linear mixed effects modeling to 
determine a value and a standard error for each of a 
plurality of allometric parameters ; 

(b) inputting said value and said standard error of 
at least one of said plurality of allometric parameters 
into a stochastic pharmacokinetic model, wherein said 
allometric parameter is designated as a random variable; 
and 

(c) using said stochastic pharmacokinetic model to 
computationally simulate administration of said compound 
to said human to produce said distribution of said 
outcome . 

[0023] "Non-linear mixed effects modeling" as 
used herein means a statistical modeling method used to 
describe, e.g., longitudinal data, wherein the model 
comprises both, fixed and random effects, and wherein the 
model is non-linear. See, e.g., "Mixed-Effects Models in 
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S and S-Plus, " Statistics and Computing, Ed. Jose C. 
Pinheiro and Douglas M. Bates, Springer-Verlag New York, 
Inc (2000), incorporated herein by reference. 

[0024] "Standard error" as used herein is a 
measure of the imprecision associated with the value of 
an allometric parameter. The standard error is outputted 
by the non-linear mixed effects modeling. 

[0025] "Allometric parameter" as used herein 
means a scaling parameter produced by correlating a 
pharmacokinetic parameter to body weight across one or 
more animal species. 

[0026] The allometric scaling equation typically 

is : 

Y = a(BW) p (1) ; 

wherein Y is a pharmacokinetic parameter such as 
clearance, volume of distribution, etc., and BW is the 
body weight of the animal studied. a and (3 are 
allometric parameters. Each allometric parameter has a 
non-zero value and a standard error associated with that 
value . 

[0027] Examples of such pharmacokinetic 
parameters are clearance, volume of distribution, and 
inter-compartmental clearance. Thus, each such 
pharmacokinetic parameter has an associated pair of 
allometric parameters, oc and (3, based on the above 
equation (1) . 

[0028] "Stochastic pharmacokinetic model" as 
used herein means a mathematical model for simulation of 
a clinical trial of a compound under investigation, 
wherein at least one of the allometric parameters is 
designated as a random variable. 
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[0029] Step (a) : 
In the first step of the present invention, 
typically, an investigational compound is administered to 
a plurality of animals. Typically, more than one species 
are used in the evaluation. One or more pharmacokinetic 
parameters are determined and the data are correlated 
with the body weight of the animal and allometrically 
scaled across the plurality of animals. 

[0030] For example, following administration of 
a compound to an animal, the compound concentration in, 
e.g., plasma, is measured at various time points post 
dose. An equation describing one or more suitable 
pharmacokinetic models (e.g. one-compartment model or a 
multi -compartment model) is fitted to these data using, 
e.g., non-linear mixed effects modeling. That equation 
is parameterized using suitable pharmacokinetic model 
parameters such as clearance and/or volume of 
distribution. These pharmacokinetic parameters are 
expressed as a function of body weight, and the 
associated allometric parameters are then determined. 

[0031] Step (b) 
In the second step of the present invention, at 
least one allometric parameter is input as a random 
variable in a stochastic pharmacokinetic model. 

[0032] The use of allometric parameters as 
random variables in a stochastic approach to human 
pharmacokinetic predictions is a key feature of the 
present invention. In the prior art, the allometric 
parameters, e.g., the relationship between a 
pharmacokinetic parameter such as clearance and body 
weight, have been treated as fixed and not subject to 
uncertainty despite the known limitations of using animal 
data to predict human data. Therefore, extrapolation of 



animal pharmacokinetics to humans is subject to varying 
amounts of uncertainty, depending on the animal species 
studied, the compound under investigation, and the study 
designs employed. However, applicant's invention uses a 
stochastic approach to interspecies allometric scaling. 
One measure of the degree of uncertainty in the 
allometric scaling is the standard errors in the 
estimates of the allometric parameters. When the 
allometric scaling relationship correlates well with the 
animal data, then the standard error is small. When the 
allometric scaling relationship correlates poorly with 
the animal data, then the standard error is large. Thus, 
applicants have recognized, for the first time, that the 
standard error in the allometric parameter estimates can 
be used to quantify the uncertainty when extrapolating 
the animal data for human pharmacokinetic /PD 
predictions. This allows a more informed selection of 
doses to investigate in the first study representing 
administration of the investigational compound to human 
subjects . 

[0033] Typically, each allometric parameter that 
is input as a random variable has an initial mean value 
equal to the value of said allometric parameter, and a 
standard deviation equal to said standard error of said 
allometric parameter. The initial mean value and the 
standard error are as determined from the allometric 
scaling of the animal data in step (a) . 

[0034] According to a preferred embodiment, each 
of the allometric parameter is designated as a random 
variable with a specified probability density in the 
stochastic pharmacokinetic model. Probability density, 
as used herein, describes the probability of observing a 
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specific value predicated on a certain probability 
distribution. 

[0035] According to another preferred 
embodiment, the specified probability density of each 
random variable independently corresponds to either a 
normal distribution or a lognormal distribution. 

[0036] Step (c) 

[0037] In the third step of the present 
invention, the stochastic pharmacokinetic model is used 
to computationally simulate the administration of the 
compound to at least one human. Typically, a plurality 
of administrations to a group of humans is simulated. 
The simulation produces the expected distribution of one 
or more outcomes of the administration of the compound to 
humans . 

[0038] Typically, such a computational 
simulation allows for variability in the distribution of 
human body weight and random inter- individual 
variability. However, applicants 1 invention allows for 
building into such a simulation the variability 
associated with the allometric scaling across the animal 
species. Thus, the distribution outcome produced by the 
simulation reflects, in part, the uncertainty inherent in 
the allometric scaling across the animal species. 
Consequently, the distribution outcome produced by the 
present invention is more useful in predicting the range 
of possible and probable pharmacokinetic /pharmacodynamic 
responses that might arise from administration of the 
compound under investigation to human subjects. 

[0039] Commercial software packages may be 
readily employed for conducting such simulations. An 
example of such a commercial software package is 
Pharsight® Trial Simulator™. See, also, e.g., "Computer- 
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assisted drug development (CADD) : an emerging technology 
for designing first- time-in-man and proof -of -concept 
studies from preclinical experiments," Gomeni , R. et al. , 
European Journal of Pharmaceutical Sciences, 1J3 (2001) 
2 61-270, incorporated herein by reference. 

[0040] Preferred outcomes produced by the 
present invention include, e.g., AUC, Cmax, C24, Cavg, 
Cmin, or a pharmacodynamic response. More preferred 
outcomes are AUC, C av g or Cmax. A preferred 
pharmacodynamic response outcome is, e.g., inhibition or 
stimulation of a biological target. 

[0041] Selection of Dose Levels 

[0042] The distribution of an outcome produced 
by the present invention is useful for selecting the 
initial (first) dose to be administered to humans. It is 
also useful for selecting doses to administer to humans 
in the context of a single dose, dose escalation study. 

[0043] A commonly used consideration in 
selecting an initial dose is the safety factor. The 
safety factor provides a measure of risks associated with 
administering a compound to a human for the first time. 
The safety factor is the ratio of a measure of exposure 
such as AUC associated with the no-adverse-effect-level 
and the anticipated exposure, e.g., AUC, in humans. A 
large safety factor suggests lesser risks in 
administering a compound to a human for the first time. 

[0044] For example, based on a distribution of 
an outcome such as AUC, an anticipated safety factor can 
be determined based on the mean AUC value predicted for 
humans or, e.g., the 95 th percentile AUC value. This 
latter safety factor, based on the 95 th percentile AUC 
value, would characterize the risk to human subjects for 
individuals with extreme values of AUC. And, it would 
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also characterize the risk arising from uncertainty in 
the allometric relationships. An allometric relationship 
characterized by a high degree of uncertainty, i.e., 
large standard errors, will result in a broad 
distribution of AUC values. Therefore, the safety factor 
determined from the 9 5 th percentile anticipated AUC would 
be much smaller resulting in a more conservative 
selection for the initial dose. 

[0045] Thus, according to an alternative 
embodiment, the present invention provides a method of 
determining a priori an optimal dose range for 
administration of a compound for the first time to 
humans . 

[0046] The optimal dose range provides a minimum 
dose level and a maximum dose level to be investigated in 
humans. The minimum dose level is an amount of the 
compound under investigation required for a desired level 
of safety and, optionally, a desired minimum 
pharmacological effect in a human. The maximum dose 
level is an amount of the compound under investigation 
required for a desired level of safety and, optionally, a 
desired maximum pharmacological effect in a human. 

[0047] Such an administration of the compound 
for the first time to humans can be part of a single dose 
study, wherein a single dose amount from within the 
optimal dose range is selected for administration to one 
or more humans. Alternatively, such an administration 
can be part of an escalation study, wherein one or more 
humans are administered a series of escalating dose 
amounts. In such an escalation study, the optimal dose 
range is useful in deciding the minimum and maximum dose 
levels . 
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[0048] According to an alternate embodiment, the 
present invention provides optimal doses of compounds 
under investigation, wherein the doses are determined 
according to the methods of the present invention. Such 
optimal doses are useful for administration to humans for 
the first time. 

[0049] In order that this invention be more 
fully understood, the following examples are set forth. 
These examples are for the purpose of illustration only, 
and are not to be construed as limiting the scope of the 
invention in any way. 

[0050] Example 1 

[0051] To identify the safety factors of an 
initial dose of an investigational compound, VX-702, to 
be administered to healthy human volunteers in a first- 
time-in-human study using an interspecies allometric 
scaling approach combined with a stochastic clinical 
trial simulation . 

[0052] In the first step, allometric parameters 
were determined as follows: 

• One compartment pharmacokinetic model was fitted to 
concentration- time data obtained from 11 mice, 2 
rats, 4 dogs, and 4 monkeys following i.v. 
administration of VX-702; 

• Population pharmacokinetic were described by a 
Nonlinear Mixed Effects Allometric Model (Fig, 1 and 
Fig. 2) ; 

• Fixed Effects: Population geometric means (]j)of 
pharmacokinetic parameters (clearance, CL and volume 
distribution, V) were described by a power function 
of body weight: 

Up = U(BW) UU ; 
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where P is CL or V (Table 1) ; 
• Random Effects: Population variability in 

pharmacokinetic model parameters was characterized 
by a log-normal distribution: 

Pi ~ JJVfop, (Op) ; 
where P± is the pharmacokinetic parameter value (CLj 
or Vi) for the i th animal and co P is the geometric 
standard deviation. 



Table 1 





Parameter 


Estimate 


Standard 
Error 


Clearance 


<*CL 


0.239 


0.0359 


Pcl 


0.603 


0.0419 




0.451 


0.118 


CD 




1.17 


0.125 


E 

O 


Pv 


0.961 


0.0609 


> 


COv 


0.402 


0.0478 



[0053] In the second step, the above allometric 
parameters, inter alia., were input into a stochastic 
pharmacokinetic model as random variables (Fig. 3) : 

• Allometric parameter estimates & standard errors 
were input into the model ; 

• Assumptions of typical bioavailability (25, 50%) and 
absorption rate constant (0.1, 0.4 hr" 1 ) were input; 

• Typical distribution of body weights in study 
population was input. 

[0054] In the third step, the administration of 
VX-7 02 to human subjects was computationally simulated 
(Fig* 3) . 250 trials of 12 human subjects were simulated 
for various starting doses: 
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• 250 trials of 12 subjects each simulated for various 
starting doses (Results presented for selected 
starting dose of 2.5 mg) ; 

• Allometric parameter values randomly selected from 
distribution described by mean and standard error ; 

• Subject body weight and intersubject variability 
randomly selected from specified distributions. 

The output of the computational simulation consisted of 
the distribution of the metrics that measure the 
exposure levels of the subjects to VX-702. 
Specifically, the output consisted of: 

• Empirical distribution of exposure metrics (AUCO-oo, 
Cmax, C 2 4 hours) (Fig. 4, Fig. 5, and Fig. 6) ; 

• The predicted safety factors for humans at each dose 
level (Table 2) ; and 

• Pharmacodynamic response: ex vivo LPS stimulated 
cytokine production . 
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Table 2 

(A) 



(B) 









2.5 


5 


10 


20 


40 


80 \ 


Scenario 


%F 


KA 


mg 


mg 


mg 


mg 


mg 


mg 


1 


50% 


0.4 


61 


29 


15 


7 


4 


2 


2 


50% 


0.1 


67 


32 


16 


8 


4 


2 


3 


25% 


0.4 


132 


61 


29 


14 


7 


4 










2.5 


5 


10 


20 


40 


80 


Scenario 


% F 


KA 


mg 


mg 


mg 


mg 


mg 


mg 


1 


50% 


0.4 


26 


12 


6 


3 


2 


1 


2 


50% 


0.1 


27 


13 


7 


3 


2 


1 


3 


25% 


0.4 


57 


26 


! 12 


6 


3 


I 2 



The predicted safety margins at each dose level are calculated as the ratio of the 
AUCo-v from the 28-day monkey toxicology study to the predicted median (A) or 95* 
percentile (B) AUCo. ¥ in humans. 
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